Majorana fermions of a two-dimensional p x +ip y superconductor 
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To investigate Majorana fermionic excitations of a p x + ip y superconductor, the Bogoliubov-de- 
Gennes equation is solved on a sphere for two cases: (i) a vortex-antivortex pair at opposite poles 
and (ii) an edge near the south pole and an antivortex at the north pole. The vortex cores support 
a state of two Majorana fermions, the energy of which decreases exponentially with the radius 
of the sphere, independently of a moderate disorder potential. The tunneling conductance of an 
electron into the superconductor near the position of a vortex is computed for finite temperature, 
and is compared to the case of an s- wave superconductor. The zero bias conductance peak of the 
antivortex is half that of the vortex. This effect can be used as a probe of the order parameter 
symmetry, and as a direct measurement of the Majorana fermion. 
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I. INTRODUCTION 

In the last few years it has appeared increasingly likely 
that nontrivial (or non-Abelian) topological phases of 
matter pQ may be produced in the laboratory. In fact, it 
is quite possible that these phases of matter have been 
produced, although current experiments still leave room 
for doubt. Interest in such phases of matter is driven to 
a large extent by their possible application for building 
naturally error resistant, so-called "topological" quantum 
computers [1 j. Among such topological phases of matter, 
perhaps the simplest is of the "Ising" or SU(2) 2 class pQ, 
which correspond to chiral p x -\-ip y BCS paired supercon- 
ductors [2]. 

There are several possible physical systems where 
Px+Wy pairing is believed to be realized, including the 
A phase of superfluid 3 He [3 ( 3 HeA), the exotic super- 
conductor Sr 2 Ru04[4 , and the v = 5/2 quantum Hall 
state [5j[6]. In addition, there have been recent propos- 
als to realize p x + ip y pairing in cold fermion gases [7] . 
For quantum information processing applications, two- 
dimensionality (or at least quasi-two-dimensionality) is 
necessary. This is certainly the case in quantum Hall 
systems, and may also be achievable for Sr2Ru04 (which 
is a layered structure), for 3 HeA films, and also poten- 
tially in cold atomic systems. For the purpose of this 
paper we will assume quasi-two-dimensionality, although 
some parts of our results are more general. 

In these (weak coupling) p x +ip y systems, certain types 
of vortices (quasiparticles in the quantum Hall context 
[8]) are believed to carry zero energy Majorana fermions 
[HI [9] , which are topologically protected degrees of free- 
dom. In Sr2Ru04 and 3 HeA the vortices that carry the 
Majorana fermions are the so-called half-quantum vor- 
tices, which can be thought of as a vortex in the order 
parameter of one spin species, without a vortex of the op- 
posite species [10,. (Note that in spin-polarized p x -\-ip y 
systems, including proposed atomic gas realizations or 
the 5/2 state, there is no half-quantum vortex and the 



full quantum vortex carries the Majorana fermion). 

A Majorana fermion is an operator which satisfies 
the fermionic anti-commutation relation {rf (x),r)(x f )} — 
2S(x — cc'), but equals to its own hermitian conjugate 
r)\ = r\. Therefore rf = rf 2 = r^rj = 1. A fermion 
occupation number state can be defined as a linear 
combination of two Majorana fermions localized in two 
distinctive vortices ip(x) = (rji(x) + ir]j(x))/y/2 and 
ip^(x) = (rji(x) — ir]j(x))/V2. This operator satisfies the 
usual fermionic relations {^(x), ip(x f )} = S(x — x') and 
-0t2 _ ^2 _ q ghgji ca ]j suc h a fermionic occupation 
number state a Majorana state. The energy of the Majo- 
rana state reflects the exponentially small hybridization 
between the wavefunctions of the two localized Majorana 
fermions. 

The zero energy of the Majorana fermion in a single 
vortex is believed to be topologically protected against 
weak disorder [81. But in the case of more than one vor- 
tex, an experimentally relevant question is whether the 
exponential localization and hybridization of the Majo- 
rana state is a property only of a clean system, or is it 
robust against the inclusion of disorder. We find that 
these properties survive even in the presence of a moder- 
ate disorder. 

Another unique property of the p x Mp y order parameter 
is the existence of low energy chiral states, which are 
localized along the edge of the sample [8]. If a single 
vortex is present, the Majorana state is split between the 
vortex core and the edge. 

Let us suppose that in one of the above systems, the 
relevant Majorana- fermion-carrying vortex has been cre- 
ated [11 j. An important next step would be to design an 
experiment to demonstrate that the Majorana fermion 
is present in such a vortex [12]. In the case of Sr2Ru04, 
one obvious experiment would be an energy-resolved tun- 
neling experiment, which measures the local density of 
states (LDOS) [T3] , An observation of a localized mode 
at precisely zero energy would be direct evidence of the 
Majorana fermion. For cold atoms, an analogous ex- 



2 



periment for observing the LDOS would be an energy- 
resolved local particle annihilation experiment. For the 
other realizations of p x + ip y order it is not as clear how 
such an experiment would be performed [T4] . 

In principle such tunneling experiments could pro- 
vide definitive evidence for the Majorana fermion. How- 
ever, in practice they may be prohibitively difficult. In 
the vortex, there will exist sub-gap bound states in 
the core known as Caroli-de-Gennes-Matricon (CdGM) 
states p~5j[l6]. The spacing between the CdGM states is 
approximately e c = A 2 /ep, where Ao is the gap (pre- 
sumably on order of the critical temperature) and ep 
is the Fermi energy. Since the experimentally observed 
tunneling spectrum will be smeared by the temperature, 
this tunneling experiment would naively only have a clear 
signature for T < e c . Unfortunately such low tempera- 
tures could potentially be unattainable in any of the pro- 
posed realizations (e c « 7fiK in 3 HeA, and < O.lmK in 
Sr 2 Ru0 4 ). 

We find that within the reachable temperature region, 
e c < T < Aq, the central peak of the smeared LDOS of 
the antivortex is half the height of the peak of the vortex. 
We shall see that this distinction is clear evidence of the 
Px + ipy symmetry of the order parameter, and of the 
existence of the Majorana fermion. 

Generally speaking, a physical asymmetry between a 
vortex and an antivortex can occur only in superconduc- 
tors which break time reversal symmetry, such as p x +ip y 
superconductors [17] . The order parameter of such super- 
conductors involves internal angular momentum, which 
is interlaced with the angular momentum of the vortex 
according to theirs relative directions [18]. 

This paper is organized as follows: 

In Sec.[n]we implement the p x Mp y superconductor on a 
sphere, with vortex-antivortex pair at the poles [19 . Us- 
ing monopole harmonics functions, we numerically solve 
the Bogoliubov de-Gennes (BdG) equation, and get the 
full BdG spectrum. 



In Sec. Ill we test the exponential decay of the Majo- 
rana state energy as a function of the distance between 
the vortices in the presence of disorder. We find it is 
unaffected, even in the presence of a moderate disorder. 



In Sec. [IV] we put an edge around the south pole of the 
sphere, and observe the edge excitations and their linear 
dispersion. 

In Sec. [V] we calculate the tunneling conductance of an 
electron into the superconductor near the position of a 
vortex or an antivortex at zero temperature and at ele- 
vated temperature. We find an asymmetry effect in the 
zero bias conductance between the vortex and the an- 
tivortex, which one can use as "smoking gun" evidence 
of the existence of the Majorana fermion. We compare it 
to the tunneling spectrum for a regular s- wave supercon- 
ductor, to support this conclusion. Our analysis shows 
that this effect will occur for any spin polarized chiral su- 
perconductor (chiral-p, chiral- d, . . . ). However, it occurs 
for the single vortex only for the chiral- p case, whereas it 
occurs for the double vortex for chiral- and correspond- 



ingly higher vortices for higher pairing symmetries. Some 
of the results of this paper were recently published in a 
short format [20l. 



II. BdG THEORY ON A SPHERE 

Consider a two-dimensional uniform p x +ip y supercon- 
ductor of spinless fermions. The excitation spectrum is 
given by the BdG equation [21] 



T - e F + W A 

At -(r-e F + W) 



En 



where T is the kinetic energy operator, and is the 
Fermi energy. W denotes electrostatic potential, and will 
be zero in this section. A is the order parameter field, 
which according to the p x +ip y symmetry is of the form 



A(x - x') = A 



(x-x') 2 



(2) 



^—[(x-x')+i(y-y')}e V . 



87ri£ p kp 

Here A is the pairing amplitude, and kp is the Fermi 
wavevector, given by ep = k F 2 /2m* , where m* is the 
electron effective mass. £ p is the pairing range, which is 
usually taken to be zero for simplicity, whereas in quan- 
tum Hall systems it is comparable to the magnetic length. 
In Fourier space 



A k = A fcx ^ fcy e- fc V, 



(3) 



which shows that Aq is approximately the energy gap. 
Although we have written down a special case of p x +ip y 
order parameter, the k — > part is universal. 

All the calculations in this paper are considered to be 
in the short range pairing limit, where kp^ p <C 1. How- 
ever, the form Eq. [3] remains acceptable up to kp^ p ~ 1, 
although due to the exponential factor the amplitude is 
highly reduced (but can be compensated by multiplying 
by a factor of e kr ^ ) . We have checked that none of our 
results change substantially even for kp^ p ~ 1. 

A vortex (+) and an antivortex (— ) are described by 
the order parameters 



A ± (x,x') = A(x-x')f v (f/Oe ±i$ , 



(4) 



where f, </> denote the polar coordinates of the pair center 
of mass (x + x')/2. The amplitude profile of the vortex 
fv(x) vanishes at the origin, and approaches unity at 
x ^> 1. £ is Pippard's coherence length [22] 



2e F 
nAokp 



(5) 



Note that for the vortex A + , the relative and the cen- 
ter of mass angular momenta are aligned, while for the 
antivortex A_, they have opposite chirality. 
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We implement the BdG equation on a sphere of ra- 
dius R, parameterized by the unit vector ft = (0,0). 
The spherical geometry has two important advantages: 
(i) It has no boundaries, which strongly aff ect the low 
energy spectrum (as will be discussed in Sec. IV). (ii) It 
enables the use of monopole harmonics functions as a ba- 
sis, which appears to be very convenient for the p x +ip y 
pairing. However, the spherical symmetry forces us to 
consider an antipodal vortex-antivortex pair. We set such 
a vortex-antivortex pair in the north and south poles, re- 
spectively, see Fig. [I] The azimuthal symmetry of this 
configuration conserves the azimuthal angular momen- 
tum, which greatly reduces the computational difficulty 
of the BdG diagonalization. 

The spherical geometry is used in this paper for cal- 
culations on finite size geometry, where the physical 
limit of far separated vortex-antivortex pair in the two- 
dimensional plane is approached in the R — > oo limit. 
Since the quantities we are interested in will be found 
to decay exponentially fast with our finite sphere 

calculations will be relevant at values of R/£ which are 
not enormously large, i.e. in moderate vortices density. 
Of course, our calculations would be even more directly 
relevant to thin spherical shells which might be experi- 
mentally created on small spherical substrates. 

The order parameter field on the sphere is taken to be 
of the following form [23] 



antivortex 



Ay («,«') 

A p (n,n') 



a 



A p (n,n')*v(n), (6) 

Aq 

(4<p 2 )(Z i , + |) 

x(af3 f -f3a f )\aa f * + pp'*\ 2 WM 2 , (7) 
cos(<9/2), 

sin(0/2)e-^, (8) 



where A p gives the pairing of the particles, and is con- 
structed by the spinor functions a and (3. The (...) factor 
in A p acquires a 2ir phase winding when fi encircles 1Y , 
which describes p x + W y pairing. The | . . . | factor keeps 
the particles within relative distance £ p . lp is the Fermi 
angular momentum, given by 

If(If + 1) 



2mR 2 



(9) 



In this way A p reproduces A (Eq.|2| in the large R limit. 

Fy(f2) describes the vorticity of the pair center of mass 
Cl = (ft + f2')/2. We choose Fy to describe an antivor- 
tex on the north pole and a vortex on the south pole, de- 
picted in Fig. [I] For the vortex pair field, we use (without 
self consistency) the approximate solution of the Gross- 
Pitaevskii equation of a vortex [24] 



sin 0-R/t 



J<t> 



^l + (sin#- R/0 2 



(10) 



with £ p < £ <C R for simplicity. 

We expand the order parameter as a series of monopole 
harmonics [25] Y q i m , of q = — | and Z,m half integers 




vortex 

FIG. 1: A vortex-antivortex pair of the 

Px ~r ^Py supercon- 
ductor on the sphere, described by Eq. (J6|. Thin black lines 
represent the current flow. Wide green arrows represent the 
pair relative angular momentum. £ p is the pairing range. £ is 
the coherence length, which determines the vortex core size. 



J. These monopole harmonics represent eigenstates 
of a particle on a sphere in a radial magnetic field with 
2q = — 1 flux quanta penetrating into the sphere (the — 
sign was chosen for correspondence with the composite 
fermion picture of the v — 5/2 state). In this basis, the 
BdG equation is represented as a matrix 



Tim, I' m' 



A v - 

lm,l'm' 



— Ern 



^n,l' rn' 
njfh 



V, 



(ii) 



with summation over primed indices, and with the fol- 
lowing matrix elements: 



Timj'm' — $ll>$mm' 



1(1 + 1) 



If(If + 1) 



x lm,l' 



-W,l-ra^0 



-(2Z + 1)(2Z' + 1) 

(r>i + (-1)'+'' A') v / 2iTT £ 



(12) 
(13) 



V L 



I V L 
-m m — 1 1 



Di 



(14) 



The matrix element A^ Vrn , (Eqs. 13fll4 ) is justified in 
Appendix [a] (A definition for f£ can also be found 
there). As expected by the azimuthal symmetry, m is 
a good quantum number. This dramatically simplifies 
the numerics, since the BdG matrix can be diagonalized 
for each m separately. 



Diagonalizing Eq. ( 11 ) for each m produces a set of en- 
ergies E n:7ri and corresponding eigenvectors w n ,im^n^m- 



4 



The resultant BdG wavefunctions on the sphere are 

U n ,m( n ) = ^2 U n,lmY-l,l, m ( n ) ( 15 ) 
I 
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FIG. 2: BdG spectrum E n ,m, of the vortex pair on the sphere, 
depicting the CdGM core states. The inset shows that their 
double degeneracies are split by weak tunneling between the 
poles. The positive energy member of the doublet which sad- 
dles zero energy at m = |, is the Majorana state shared by 
the vortex and antivortex. 

In Fig. [2] we depict the BdG spectrum of the vortex 
pair state as a function of m. The spectrum shows an 
expected symmetry of the BdG equation, which implies 
that for every eigenvector (u ni v n ) with energy E n , the 
vector is also an eigenvector with energy — E n . 

Hence according to Eqs. (15)-([l6|), 



On' 



-ra+1 



(17) 



En,', 



ra+1- 



The continuum states above the gap \E n \ > Ao are 
extended, while the branch that approaches zero is the 
p x -HPy version of the CdGM core states. Their dispersion 
is 



El 



(18) 



and their number is of order ep/^o [16] . 

As seen in the inset of Fig. [2j each CdGM state is al- 
most doubly degenerate. The splitting represents weak 
tunneling between the north and south pole cores, and 
decreases exponentially with the radius of the sphere 
SE^ ~ e~ R ^ for R ^> £. In particular, the lowest 
positive energy E +, and its negative companion E - 
(both at m = approach zero as e -jR ^, is shown in 

2 = \v A -m 2 

Fig. [4] The 

wavefunctions are symmetric and anti-symmetric super- 
positions of the north and south localized core states. 



Fig. [3] The probability densities \uq+(9)\ z 
and \u -(0)\ 2 = \v +(0)\ 2 are depicted in 
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FIG. 3: Energy of Majorana state E + as a function of sphere 
radius R. The exponentially decreasing energy indicates tun- 
nel splitting between vortex and antivortex core states. 




FIG. 4: Probability densities of Majorana states \uo\ 2 ver- 
sus latitude on the sphere 0. The smooth parts of the Ma- 
jorana state wavefunctions u +(0) and u -(0) are approxi- 
mately symmetric and antisymmetric with respect to reflec- 
tion about the equator = tt/2. Both \u + | 2 (blue) and \u - | 2 
(green) show the exponential localization in the (anti)vortex 
cores. 



In the infinite sphere limit E { 



0+ 



0, and 



the wavefunctions uq(Q) ~ ^o(^) are equally split be- 
tween the north and south poles. The corresponding 
BdG quasiparticle field operator near each pole r](fl) = 
uo(ft)^(ft) + vo(Q)^ (ft) , is a Majorana fermion oper- 
ator 77 « rf. 

The asymptotic behavior of the wavefunctions uo(r) in 
the plane are [TJ [T6] 



uq(x) 



J (k F r)e-^f rd ( r '/^ f ^ r '/^ antivortex 
Ji (k F r)e- i f r d ( r ' / fAr / e ^ vor tex 

_ (19) 
where / v was defined in Eq. Q. We confirmed this 
asymptotic behavior, as seen in Fig. |5| 



The physical reason behind the difference in Eq. (19) 



between the vortex and the antivortex is that the phase 
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antivortex 




FIG. 5: Probability distributions of Majorana state |^ +| 2 
as a function of the distance from the (anti) vortex center r, 
at the cores vicinities. The wavef unctions have exponentially 
small weights near the equator. The asymptotic (red) and the 
numeric (blue) wavefunctions show excellent agreement. 



winding of the order parameter is determined by the sum 
of the vorticity and the relative angular momentum. For 
a vortex the vorticity and the angular momentum are 
aligned, which yields the phase winding e 2 *^. In that 
case the condition for Majorana fermion solution of the 
BdG equation, u(x) = v*(x), can be fulfilled only by 
uq(x) = uo(r)e l( P , where uo(r) vanishes at the origin (the 
J\ (r) behavior is expected, due to the azimuthal angular 
momentum 1). By contrast, for an antivortex the order 
parameter is real, since the vorticity and the relative an- 
gular momentum cancel each other. Thus the Majorana 
fermion solution is real and radial uo(x) = uq(v), and can 
be finite at the origin. (Here, Jo(r) is expected, since the 
angular momentum is zero). We will see in Sec. |v| that 
this difference is crucial for an experimental signature of 
the Majorana state. 

The excellent agreement between the asymptotic and 
numerical wavefunctions, as seen in Fig. |5j reveals the 
underlying physics of the core states. The radial profile 
of the order parameter serves only as a confinement, and 
therefore determines only the exponential decay of the 
wavefunction (and very weakly the spacing between the 
core states). The short range part of the wavefunction 
is controlled by a Bessel function, which is determined 
only by the symmetry and the total angular momentum. 



Furthermore, the specific pairing functions and pairing 
range do not play any apparent role. 

s-wave . In order to emphasize the unique behavior 
of the chiral p-wave order parameter, we compare it to a 
regular s-wave order parameter, characterized with the 
same set of physical parameters. 

An order parameter with s- wave symmetry in a plane, 
in the presence of a vortex, will have the form 



A s ± (x,x') = A s (x — x')f v (f/£)e 



A s (x-x') = A 



(20) 
(21) 



where / v is the same as in Eq. A s provides the 

s-wave pairing, with the same pairing range £ p of the 
Px + Wy order parameter. 

Implementing the order parameter on a sphere with a 
vortex-antivortex pair at the poles is quite similar to the 
Px + Wy case: 



A^(n.n') = A^(n,n')F v (Ci), 



(22) 



Aq 

4tiV 



\aa'* + l3p'*\ 2(R >^\ (23) 



where Fy is unchanged. 

Without the chirality, the order parameter can be ex- 
panded in the regular spherical harmonics Yi m basis, with 
I, m integers. In this basis, the matrix forms of the kinetic 
term and the order parameter are 



Tl m ,l'm' — 8u'5 mTn i €p 



1(1 + 1) 
l F (l F + 1) 



1 , 



(24) 



lm,l' m' 



= -<y rn /,i_ rn A J I fc(2Z + l)(2Z' + l) (25) 



I V L 




I V L 
-m m — 1 1 



D 



-i 2 (i P /R) 2 



(26) 



Our method of computing Vrn , is almost the same as 
that for Aim^m' (Appendix [AJj , but with q set to and 
integer Z,ra. Df is given essentially by Eq. ( |A7[) , and 
is approximated in a way analogous to D\ (Eq. |14|), but 
without the chirality factor l/lp. 

The BdG spectrum of the s-wave system is generally 
similar to that of the p x +ip y system, as seen in Fig. [6j 
which was plotted with the same parameters that were 
used to plot Fig. [5] Both the extended states, with \E n \ > 
Aq, and the CdGM states, with « (| — m)e c , are 
present. Here m is an integer, the energies of the CdGM 
states are shifted by \e c compared to the p x + ip y case, 
and the Majorana state is absent. 

Far from the vortex core, the s-wave and the p x +ip y 
show almost the same gapped spectrum. Fourier trans- 
forming the BdG equation in the plane (Eqs. filmland 21 
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FIG. 6: The BdG spectrum E n , m , of an s-wave supercon- 
ductor, with the same physical parameters of the p x +ip y in 
Fig. [2] The distinction from the p x +ip y case can be noticed in 
the inset, where the energies of the CdGM states are shifted 
by half the level spacing, therefore lacking a Majorana state 
which saddles zero energy. 



see Eq. [3| , and diagonalizing the BdG matrix, gives 



,Px~HPy 

k 



Ei 



(k 2 /2m* - e F ) 2 + (A ^e- fe2 ^ 2 ) 2 ,(27) 
kp 

(k 2 /2m* - e F ) 2 + (A e- fe2 ^ 2 ) 2 , (28) 



where k = {k 2 x + k 2 ) 1/2 . 

The small difference between the two cannot be dis- 
tinguished in a tunneling experiment. Nevertheless, the 
disparity in the CdGM states, and especially the exis- 
tence of the Majorana state, may be observed by taking 
a difference in the tunneling LDOS, as we will see in 
Sec. El 

The last point to be mentioned is the insensitivity of 
the above results to the choice of hamiltonian parameters. 
According to Eqs. (12)-(14) there are three free dimen- 



sionless physical parameters: Ao/e^, lp and while 
£/R is determined by lp • Ao/e^ (according to Eq. 5]). 
However, the asymptotic wavefunction (Eq. [l9| ) implies 
that the exponential envelope is controlled by £. More- 
over, the energy of the Majorana state is also determined 
by £• We confirmed this numerically, for the physical 
regime Aq <C ep. hp, which is approximately IpR, con- 
trols the oscillation frequency of the wavefunction, as 
seen in Eq. [19] The gap A and cp determine the en- 
ergy spacing of the core states e c and their number. £ p 
has negligible effect, up to kp ■ ^ p < 27r, meaning that 
the pairing range has no significant effect either on the 
wavefunctions or on the energies of the core states. Ad- 
ditionally, modifying the vortex profile of Eq. (10) - to, 
for example, tanh(sin# • R/£) - only modifies the wave- 



Ill. DISORDER 

Disorder mixes states with different angular momen- 
tum m, thus making the BdG equation extremely hard 
to solve analytically. We study numerically the effect 
of disorder on the Majorana state, by the addition of a 
white noise random real potential to the BdG equation 
(Eq.0, given by 



l,m 



(29) 



where the i^ m 's are independently identically dis- 
tributed, with wi-m = (— l) m iu z * m . I \ is an ultra violet 
cut-off. We want the potential to be independent of the 
radius of the sphere via /a, therefore we take the real and 
imaginary parts of wi m to be uniformly distributed in the 

, where Wq is given in units of 



interval r 

La 

energy. This gives 



6ttWo y 67tWq 



Ia 



(W 2 (ft)) = w 



(30) 



as shown in Appendix [B] and demonstrated in Fig. [7] 

The disorder breaks the azimuthal symmetry, so that 
m is no longer a good quantum number. The matrix 
elements of W(f2) are 



Wi 



ml' m' 



(-i) m -V^( 2z+i )( 2Z,+i ) 

^ L-rri' —rn 

(31) 

L 

I V L 



I V L 

1 - 1 

2 2 u 



functions slightly inside the cores, according to Eq. 19 



as derived in Appendix |B| 

Fig. [8] depicts the disorder averaged energy of the Ma- 
jorana state E + versus R for increasing Wo- It can be 
seen that both the average energy (solid) and the stan- 
dard deviation (error bars) decay exponentially in the 
regime Wo < ^f- Thus we conclude that the exponen- 
tial drop of the Majorana state energy with increasing 
system size survives moderate disorder. 



IV. SYSTEM WITH AN EDGE 

The p x +ipy state has broken time-reversal symmetry, 
implying there are chiral modes which are exponentially 
localized at the edge of the sample. For a disk of radius 
R, the energies of the edge states are expected to have 
the form E^ ge oc m/R, where the angular momentum 
m is a half integer, due to the anti-periodic boundary 
condition of the spin polarized fermions [3 [26] . 

In the presence of a half-quantum vortex, the boundary 
condition on the BdG wavefunctions is periodic, and m 
is an integer, with the Majorana state having the m = 
quantum number. A fermion occupation number state is 
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FIG. 7: Three white noise potentials VK(fi) on the sphere, 
taken from the same distribution of harmonics components 
wim, but of increasing high angular momentum cut-off I a 
(from top to bottom). Each time I a is multiplied by 2, while 
the disorder strength Wo remains constant. 



FIG. 8: Disorder averaged energy of Majorana state E + 
versus R/£ for increasing disorder strength Wo (from top to 
bottom). The exponential decay survives in the regime Wo < 

6F- 



created from a combination of Majorana state on the edge 
and in the vortex core. This was recently investigated 
numerically [27] . 

In order to create an edge at latitude Be, we add a 
strong potential of the form 



W E {0) 



2e F 



e (0 E -e)2R/£ + i 
= 5>f>Lo(ft), 



(32) 



which defines wf. Note that at the edge We(0e) = cp, 
which sets the density to zero. The width of the potential 
was chosen to be the longer length scale £. 

For a uniform superconductor with an edge, the order 



parameter is of the form 



Atf(n,n') 

F E {0) 



A p (W)F E (Cl), (33) 

tanh {{0 E - 0)R/£) 0<B <B E 
6 E < < 7T 







£/^o(n), 



(34) 



where the pairing A p (n,n f ) is defined in Eq. Q, and 
the envelope 2^(fi) equals to zero for Be < < 7T due 
to self consistency. The matrix elements of We and A e 
appear in Appendix [Cj which also shows that m is a good 
quantum number. 

Fig. [9] depicts the BdG spectrum of an edge without 
a vortex. The states above the gap \E n \ > Ao are ex- 
tended, while the branch is composed of the chiral edge 
state. However, since there is only a single edge, these 
states are not degenerate. 

In a half-infinite plane the dispersion of the edge states 
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FIG. 9: BdG spectrum E n , m , of an edge without a vortex, 
depicting edge states. The solid black line is mAo/If. Here 
R/£ = 7.6, while the fit improves for larger R/£. 



Note that the vortex envelope Fy E is no longer symmetric 
with respect to the equat or, an d so includes f^ E / for 
even L's, unlike Fy (Eq. |A14 ). The matrix elements of 
Aye appear in Appendix|C| 




«^^>^^ ' Majorana 

< X < < > * : 2 v-vi v »* : i . !■ . t I 

I** t* * * ♦ *♦ r** t *t ** £ *~ ♦* **t*L^ ******* 
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20 



111 



is given by [26] E^ dge ^ Ao^, with the velocity i^d^e = 
Conversion to spherical geometry with R ^> £, gives 



m 
1^' 



(35) 



The black line in Fig. [9] depicts this approximation for 
R/£ = 7.6. The approximation improves for larger values 
of Moreover, Fig. [To] shows that ££^ e ex 1/R, as 

expected. 
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FIG. 10: Energies of the two lowest edge states without a 
vortex Em 9e , for m = | and m = |, versus £/R, show the 
decay. 



When considering the antivortex in the north pole, the 
potential remains the same, while the order parameter 
becomes 

A ra («,«') = A p (n,n f )F VE (Cl), (36) 

tanh (sin {^j^j f ) < < E 
6 E < < 7T 

(37) 



E/rni(n). 



FIG. 11: BdG spectrum E njm , of an antivortex with an 
edge, showing both the edge states and CdGM states. The 
Majorana state is present, with an exponentially small energy, 
as seen in the inset. 



The 
Fig. 



resultant BdG spectrum, E nrn , is shown in 
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The edge states and CdGM states appear in 
two separate branches. The only quasi degenerate state 
is the Majorana state at m = |. 

Fig.[l2]depicts the wavefunctions of the Majorana state 
%+(#), the first excited edge state u ed ^{&) at m = |, 
and the first excited CdGM state u c _ 1 ^ 2 (6) at m = — \. 
It can be seen that the Majorana state has almost equal 
support on the edge and in the vortex core (and is ex- 
ponentially localized in both), while the edge state and 
the CdGM state are concentrated either at the edge or 
at the vortex core, respectively. Fig. fl3] confirms that as 
expected, the energy of the Majorana state E + decays 
to zero as e _jR//< ^, while the energy of the first excited edge 
state E^j ge scales as 1/R. 

In summary, the spherical geometry enables an easy 
visualization of edge effects, by the addition of a confining 
potential. The expected edge states appear in any case, 
while the Majorana state appears only in the presence of 
a vortex. 



V. TUNNELING LOCAL DENSITY OF STATES 

The energy gap and the coherence length of a super- 
conductor can be detected by tunneling of electrons to 
its surface. For high spatial resolution, the tunneling 
also detects the low energy excitation spectrum inside a 
vortex core. In this section we show that tunneling exper- 
iment also provides direct signatures of the symmetry of 
the order parameter, and the existence of the Majorana 
state. 

At zero temperature the tunneling local density of 
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£f J -10 



FIG. 12: \og\u(0)\ 2 of the Majorana state u + (top), the 
first excited edge state (middle), and the first excited 

CdGM state u c _ 1 ^ 2 (bottom). The Majorana state is split 
between the edge and the vortex core, while in each branch 
of excitation the wavefunction is concentrated either in edge 
or in the vortex core. 



states (LDOS) is defined as [13] 
T(E, r) = V \u n (r)\ 2 S(E - E n ) + \v n (r)\ 2 S(E + E n ) 



(38) 

where r is the distance from the vortex (or antivortex) 
center. Fig. [l4| shows the LDOS at the core, for displace- 
ments r < 0.3£ and energies \E\ < 0.2Aq. The Px+Wy 
state shows a distinction between the antivortex and the 
vortex, while in the s- wave they are the same. 

The Majorana state can be easily discerned as zero en- 
ergy peaks near the centers of the vortex and antivortex 
cores of the p x +ip y superconductor. The other CdGM 
core states also appear as oscillatory peaks, with energy 
spacing e c . The difference between the vortex and an- 
tivortex of the p x + ip y LDOS is apparent: according 
to Eq. (19) the antivortex Majorana state wavefunction 
behaves like Jo, which is peaked at r = 0. By con- 
trast, the vortex Majorana state wavefunction behaves 
like Ji, which is zero at r = 0, and has a smaller peak at 
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FIG. 13: Energies of Majorana state E + (bottom) and first 
excited edge state E^ e (top) of an antivortex with an edge, 
as a function of the radius of sphere R. The energy of the 
Majorana state decays exponentially, while the energy of the 
edge state scales as 1/R. 



r « Af/4, where Xf = 27t//cf is the Fermi wavelength. 

The origin of this distinction was discussed in Sec. [TTJ 
According to that argument, in the antivortex only the 
wavefunction of the Majorana state, with m = |, is real 
and peaked at the origin. All the excitations with m ^ \ 
have a phase winding, and are therefore equal to zero 
at the center of the vortex, and have much lower peaks. 
On the other hand, in the vortex the Majorana state 
wavefunction has a phase winding of e^, and must vanish 
at the origin. But the u m (x) part of the first excited 
state, with m = — \ and Ey 2 = e c , is real and finite at 
the center of the vortex. Similarly the v m (x) part of its 
negative companion, with m = | and E^ 2 = — e c , is also 
real and finite at the center of the vortex. Therefore in 
the LDOS the first two excitations have relatively large 
peaks at r = 0, while all the other CdGM states are much 



lower, as seen in Fig. 14 



The same argument holds for the s- wave superconduc- 
tor. The phase winding of the order parameter with a 
vortex is e*^, while with an antivortex it is e~ %( ^. This sign 
of the phase winding is only a matter of convention; there 
is no distinction in the s-wave superconductor between 
the vortex and the antivortex. Physically speaking, the 
pairing of the particles does not involve any internal an- 
gular momentum, so the sign of the vorticity is meaning- 
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antivortex 




FIG. 14: Zero temperature LDOS T(E,r), in the cores of 
the p x -\-ip y antivortex (top) and vortex (middle), and of the 
s- wave vortex (bottom). e c is the spacing between CdGM 
states, and Xf is the Fermi wavelength. The peaks belong to 
the CdGM states. Notice that the zero energy Majorana state 
is maximized at the origin in the p x -\-ip y antivortex, while it 
is removed from the origin in the vortex, and is absent in the 
5-wave . 



ical potential and temperature T. In the BCS weak 
coupling regime, Ao <C e^, and therefore e c could be 
a very small temperature scale. At moderate tempera- 
tures e c < T < Ao, the peaks of Fig. [l4]are smeared on 
the energy axis (but not on the r axis). 



antivortex 




less. For such an order parameter there is no Majorana 
fermion solution, since it requires ^o(^) = iio(V)e Z( ^ 2 , 
which is not single valued. However, for the lowest state, 



with m = and En 



is real for m = 1 (and E\ 



2 c 



l (x) is real. Similarly v m {x) 
— \e c ). Therefore, the two 



lowest states are peaked at the origin, while all the higher 
excitations are equal to zero at the origin, and therefore 



much lower, in agreement with Fig. 14 



In a tunneling spectroscopy experiment (e.g. Ref. [28]), 
the tunneling conductance is measured. The conductance 
reflects the smearing of the LDOS by temperature broad- 



ening 



-,(E,r) J dE' ^^^ TiE^r), (39) 
where f{E) is the Fermi-Dirac distribution at zero chem- 



dl 

dV [ 



FIG. 15: Tunneling conductance ^(E, r), in arbitrary units, 
of the Px+ipy antivortex (top) and vortex (middle), and of 
the s- wave vortex (bottom). Ao is the energy gap, and £ is 
the coherence length, which specifies the radiuses of the cores. 
The temperature T = 0.15Ao, is about 10 times larger than 
the CdGM level spacing. The conductance of the p x + ip y 
vortex is almost identical to that of the s-wave , while the 
central peak of the antivortex is twice lower. 



The tunneling conductance at T = 0.15A = 7.5 e c , 
for displacements r < £ and energies \E\ < 1.5Aq, is 



depicted in Fig. 15 for the p x +ip y antivortex and vortex, 
and for the s-wave vortex. The conductance shows a 
central peak at E = 0, r = 0, with low broad ridges 
dispersing away to larger E,r. We see that while the 
conductance of the p x + ip y vortex is almost identical to 
the conductance of the s- wave vortex, the central peak of 



11 



the antivortex is half the height of that of the vortex. 
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FIG. 16: Zero bias conductance peak at the (anti) vortex 
core as a function of temperature T, in log-log scale. For 
e c < T < Ao the p x + ip y vortex (green) and 5- wave (blue) 
peaks are twice the height of the antivortex peak (red). 



The explanation for this factor of 2 comes from the 
origin of the zero bias peak of the conductance. Consider 
the conductance of the s-wave vortex. Since this lacks 
a zero energy CdGM state, the zero bias conductance 
vanishes at T <C e c . When the temperature is raised 
to ^e c , the two lowest CdGM states, which are peaked 



at the origin for E w ±<j>e c , broaden to create a zero 
bias peak. As the temperature rises further, this peak 
diminishes due to the broadening, as shown by the blue 



curve of Fig. 16 



The p x +ipy vortex behaves similarly to the s- wave case, 
as shown by the green curve of Fig. [l6j For T <C e c the 
zero bias conductance peak vanishes, since the Major ana 
state wavefunction has a node at the origin. By contrast, 
for T > e c a zero bias peak is created by the two broad- 
ened first excited states, which are peaked at the origin 
with E « d=e c . On the other hand, in the antivortex the 
Majorana state wavefunction is peaked exactly at the ori- 
gin. Therefore its conductance is maximal at T = 0, and 
diminishes with rising temperature, as depicted by the 
red curve of Fig. [l6j Furthermore, since in the antivor- 
tex only the Majorana state is peaked at the origin, the 
height of the zero bias peak of the antivortex is deter- 
mined by the broadening of a single state, while in the 
vortex and in the s- wave it is a sum of the broadening of 
two states. This is why the central peak in the antivor- 
tex is half the height of that found in the vortex, when 
e c < T < Ao. Thus, the asymmetry between the vortex 
and the antivortex is not only a clear fingerprint of the 
Px + Wy symmetry of the order parameter, but it is also 
a smoking gun evidence of the existence of the Majorana 
state itself. 

The distinction between the vortex and the antivortex 
spectra is based on the distinction between the wavefunc- 
tions, which may be seen on scale of the Fermi wavelength 
Xf- Therefore, as long as the spatial resolution in the 
tunneling conductance Sr is better than A F , our effect is 




FIG. 17: Effect of spatial resolution. Ratio of vortex to an- 



tivortex conductance peak heights (shown in Fig. 16 for per- 
fect resolution), for spatial resolution Sr. Xf and £ are Fermi 
wavelength and coherence lengths respectively. Temperature 
is 7.5 e c . 



observable. But when ^(E,r) (Eq. 39) is smeared over 
a length scale (Sr) 2 > X F , the ratio between the vortex 
and antivortex peak heights rapidly approaches unity as 
Sr > Xf, as shown in Fig. |TT| The requirement of such 
a high resolution, in spite of being restrictive, is feasi- 
ble in present technology of Sr ~ lA, since in Sr2Ru04 
A F ^8.3A[30]. 



Moreover, although we found in Sec. [TTT] that the Ma- 
jorana state survives moderate disorder, the suggested 
results hold provided the disorder is weak on the scale of 
Af, which is very reasonable as that limit corresponds to 
minimum conductivity and Anderson localization. 

The above analysis is based on the two-dimensionality 
of the sample. In real three-dimensional samples the con- 
ductance peaks are strongly suppressed by bulk states 
and surface imperfections [28, 31]. Nevertheless, it is the 
difference between a vortex peak and an antivortex peak 
which is sensitive to the existence of the Majorana state. 
Moreover, according to the asymptotic analysis [16 , mo- 
mentum along the vortex line (additional e %kzZ in the 
wavefunction), leaves the dispersion of the core states 
(Eq. 18) unchanged, and only modifies weakly the oscil- 



lations in the radial part of the wavefunction (through 
the oscillatory argument, which becomes \/k F 2 + k 2 r). 

Performing such an experiment is certain to be a 
challenge, although likely remains possible. One might 
imagine leaving the tip of the STM in the same position 
while reversing a weak magnetic field where the field is 
weak enough so as to flip the direction of the vortex, 
but not to overturn the chiral order parameter. Such 
an experiment would rely on an assumption that the 
vortex prefers to sit at one particular position in the 
sample — presumably due to some inhomogeneity 
or disorder in the sample — and that this preferred 
position does not change when the weak magnetic 
field is reversed. In practice, however, one would 
have to check this assumption by spatially scanning 
the STM as even a small change in the vortex posi- 
tion would have a large effect on the measured tunneling. 
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General Cases. A difference between vortex and an- 
tivortex conductance peak is expected for any chiral sym- 
metry breaking (CSB) superconductor. The important 
questions are whether this difference is observable for 
T > e c , and whether it is sensitive to the existence of 
Majorana fermions. 

For a CSB superconductor with relative angular mo- 
mentum M = 1, 2, . . . (chiral-p, chiral- . . .), in the pres- 
ence of a vortex with vorticity N = ±1, ±2, . . ., the phase 
winding of the order parameter is e ^ M+Ar ^. A Majorana 



fermion solution uq 



requires m = (M + N)/2, 



and would be single valued only if m is an integer, i.e. 
M + N is even. However, only if M + N = can the 
wavefunction of the Majorana state be finite at the center 
of the vortex, i.e. an antivortex with vorticity N = —M. 
Therefore, although a Majorana state is expected for ev- 
ery even M + N [9] , the asymmetry in the zero bias con- 
ductance peak for T > e c is expected only in the subset 
\N\ = M. Notice, that the p x +ip y vortex is the only case 
where the asymmetry occurs for a unit vorticity. 



VI. SUMMARY 

Implementing the BdG equation of a p x +ip y supercon- 
ductor on a sphere had several advantages. It allowed us 
to study the system with and without vortices and edges, 
without artifacts of finite size boundaries. We found the 
exponential decay of the Majorana state energy with re- 
spect to the intervortex distance to be robust with re- 
spect to the addition of moderate disorder strength, up to 
the Fermi energy scale. The tunneling conductance was 
plotted, and compared to the s-wave case. It was found 
that in the moderate temperature regime e c < T < Ao, 
the zero bias conductance in the center of the p x + ip y 
antivortex is half the height of that of the vortex. The 
asymmetry in the conductance peak heights is a smoking 
gun signature of the spin polarized p x +ip y superconduc- 
tor, and a direct measurement of the existence of the 
Majorana state. This idea can be generalized to order 
parameters of higher angular momentum. 



ACKNOWLEDGEMENTS 



We thank Ady Stern for useful discussions. Support 
from US - Israel Binational Science foundation and Is- 
rael Science foundation is acknowledged. AA acknowl- 
edges Aspen Center for Physics for its hospitality. HAF 
acknowledges the support of the NSF through Grant No. 
DMR-0704033. AA and SHS acknowledge the hospitality 
of the KITP, where this collaboration was initiated. 



APPENDIX A: THE MATRIX ELEMENT A^ Vrn , 



In this appendix we derive A^ Vrn , (Eqs. 13fll4 ) from 
Ay (Eqs. ^P ). 

The spinor functions a and (3 (Eq. [8| obey 

a(3 , -f3a' = 2tt (y_i i i (fl)Y_ 1 1_ i (ft') 

-y_i i_i(«)y_i i i(ft / )) , (Ai) 

22 2 V/ 2 2 2 V / 



and 



,/* +/ ^/*|2(iV£ P ) 2 



1 + cos cos 0' + sin sin 0' cos(0 - </>') ^ (jR/U 



cos - cos h sin - sin — e 

2 2 2 2 



1 + cos6 > 
—2-) 



(A2) 



where is the angle between ft and ft' . Note that for 



i + coseY" /sp/ (fl/gp)* ft2 

— ^— J «e" « 



(A3) 



which is equivalent to Gaussian pairing in the plane. 
If we denote M = [(R/£, p ) 2 ] , then 



\aa'* +/3(3'*\ 2 ( R ^ 2 



^(l+cos9) M (A4) 



M 



M 



n=0 v 7 

Any monomial x n can be expressed as a series of Legen- 
dre polynomials Pi{x) [32] 



x 



E 



(21 + l)n! 



*=„,„- 2 ,...2^(^)!(Z + n+l)!! 
which satisfy [33] 

47T 



Pl(x), (A5) 



Pz(cos0) = 27Ti i- l ) mY i-m (n)y im (no, (A6) 



m— — l 



where Yj m (IWi) are the spherical har monics. Substi- 
tuting Eqs. (A5) and (A6) into Eq. (A4), we obtain 



(A7) 



Ml 



M I 



E (-ir^- m (n)y lTO (n'), 



Z=0 m=—l 
M 



A, 



M 



E 



2 V f H_i ] i ( M - n )! (I + n + 1)!! 



n=l,l+2,... L 
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Multiplying two harmonics yields [25] 

W«)*Vl'm'(«) = ^y+l'-l-l'-m-m' 
l+V 



Alternatively, substituting k x + ik y — > (Z+ in the 
planar Ak (Eq. [3]), yields an excellent approximation for 
(A8) A 



^(2Z + 1)(2Z' + 1) (-l) ^ "^/2l 77 TI 
i"=\i-i'\ 



D, 



l+ \ e -(/ + I)2(| p/i J)2 



(A12) 



X 



Z V 



I" 



m m —m — m 



z z ; z" 

(7 g' -g - q' 



" q-\-q' ,1" ,771+771/ 

In particular 

Y olm (Sl)Y_ ii ± i(fi) = (-ir- 1 



which is valid for > 5£ p . For Z, Z^ ^> 1 we obtain 
Eq. 

For short range pairing £ p <C i?, the vorticity of the 
center of mass can be approximated by 



^(2Z + 1) ^ (-l) sT V2/' + l±2s 



Fy(O) W §(Fy(n) +Fy(n')). 

is expanded in spherical harmonics, 



(A13) 



z 



Z + s 



m ±| — m =F I 
x^-iz+ s ,m±i(") 



Z i Z + s 
u 2 2 



Z±ra + 1 
2Z + 1 



; (n)± 



F v (fl)= 53 /JIl!^), (A14) 

L=l,3,5,... 

which defines 

Substituting Eqs. ( |A1Q[ ) and ( |A14[ ) into Eq. ([6|, and 
using Eq. (|A8|), yield 



^Y 1 ^-i, ra±i (n)),(A9) 



where the last expression comes from writing the 3j sym- 
bols explicitly [33] . 

Substituting Eqs. (|Al}, ([A7|) and (|A9| into Eq. Q 
yields 



A v (n, n') = =§ ^(-i) i+i ' +L A • /j 



lm 
I'L 



x V / T i ¥ (2l + l)(2l' + l)(2L + l)(^ J ? 

; f m l l 1 )y_i Jm (n)y_i,,_ TO+1 (n') 



-m 



M 

c 

£ (-ir +1 r_ ^ +l ,_ m _i(fi)r_ 



A (n n') - Aq M 2 -M! ^ m _ M 
/\ p {ilM)- R22M+1{1 ^ + ^Ai A i+i 



m 1 _ m _ 1 ^ ^-^ / m+l(^)^-iZ-m(^ / ) 

= i E A *V m ^- i«m(n)y- i«w(n')- (Ai5) 



Iml'm' 



Ao 
i? 2 



The last line defines Aj^ z , m ,. Finally, by using 



(A10) 



r 1 v 

m" m 1 
(_l)W+«" 



(A16) 



with 



z /' 



i" 



M 2 ■ M\ 

2 M +^(l F + \) 

M-l- 

E 



Bf 1 



(AH) 



we obtain 



lm,l'm' 



1 



n=0,2,... xM-Z+i-n 



2Z + 2 + n 



J m/ ,_ m+1 A ^ T ^(2Z + 1)(2Z / + 1) (A17) 

x si (a - (-1) W ' +L A<) 



22 



(|)!(M-Z-i-n)!(2Z + n)!! 



Z Z' L 

2 2 u 



Z Z' L 
-rn m — 1 1 



where 
only for Z 



M _ 



2 K Af - i +5)(|(M-Z+|))!(M + Z + |)!! 
= M (mod 2). 



- 1 Since L is odd, we obtain Eq. (13). Note that 

A ?ki' ,-m+i = - A r,-7n+M,7n> as expected by the an- 
tisymmetry. 
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APPENDIX B: THE DISORDER POTENTIAL 

The white noise potential is defined as a series of inde- 
pendently identically distributed complex elements i^ m , 
with both real and imaginary parts uniformly distributed 



in the interval 
fore 

( w lm W l'n 



s/6ttWq V&nWo 
Ia ' *a 



Sll'&mm'dwiml 2 } 

2{\Re{w lm )\ 2 ) 



In real space 



47T 2 



Ia 



for every I < l\. There- 
(Bl) 

VJ < l A . 



^(wt m w Vm >)Y^(n)Y Vm >(n) 



l,m 
V ,m' 



In particular for the north pole 



Ia 



(B2) 



(w 2 (o = o)) = w 2 SE W« = o)| 



i=o 



A 1=0 



(B3) 



In order to calculate the matrix element Wimj'm', we 
use the identities l25l 



and 



Y q * lm (Cl) = (-iy +m Y- ql - m (Cl), (B4) 
J dnY qlm {il)Y qlVml {il)Y ql , v , m ,,{n) = 

(-i)' +i ' + V^( 2i + i )( 2i '+ i )( 2 '"+ i ) 



9 9' 9" 



z r v- 



(B5) 



Hence 



J &nY\ lm {n)w{n)Y_ ¥ , m ,{si) (B6) 
( _i r - J £ Wl „ m „ 



V'm" 



x | doyx,_ m (o)y_ |rm ,(fi)y ^(«), 



which yields Eq. ([3TJ. Note that W Vrn ^im = W^ >z , m ,, as 
expected from a real potential. 
APPENDIX C: MATRIX ELEMENTS OF EDGE 
OPERATORS 



Both the edge potential operator We (Eq. [32]) and or- 
der parameter A# (Eq. [33] ) have the following matrix 
elements: 



Imd'm' 



(5 m / > - m (-l) m -i v /^ F (2Z + l)(2^ + l) 

X ^ 2L + 1 W L ( C1 ) 



Z Z' L 

2 2 U 



Z Z' L 

m —m 



Hm,U 



rn' ,—m 



A ^(2Z + 1)(2Z' + 1) (C2) 



Z Z' L 

2 2 u 



Z Z' L 

m —rn 



Vrn , is obtained by substituting m = m' in Eq. (31) 
(due to Eq. [321 . A^ 



imi'm' i s essentially Eq. (|A17[). But 



since F E (Eq. |34| ) replaces Fy (Eq. |A14 ), m is coupled 
to —m, and /f 7 ^ for both odd and even L's. 



For the case of an edge with a vortex, the matrix ele- 
ments of the order parameter Ave (Eq. 36) are 



A VE 



Imd'm' 



5 m / > _ m+ iAo V / I ^ F (2Z + 1)(2Z' + 1) (C3) 
x £ v 7 ^! f v L E (A - (-1) W/+L A') 



Z Z ; L 

2 2 u 



Z Z' L 
-m rn — 1 1 



which are exactly Eq. ( |A17j ), with 7^ for both odd 
and even L's. 
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